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It is assumed that the quantum state that may describe a 
macroscopic system at a given instant of time is one of the 
eigenstates of the reduced density matrix calculated from the 
wave function of the system plus its environment. This im- 
plies that the above quantum state is a member of a special 
orthonormed set of states. Using a suitable Monte-Carlo sim- 
ulation, this property is shown to be consistent with the ex- 
tremely small standard deviation for the coordinates and the 
momenta of macroscopic systems. Consequences for statisti- 
cal mechanics and possible observable effects are discussed. 
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Despite of the well known conceptual difficulties , 0, 
j| , there is no compelling reason to believe that quantum 
mechanics fails when applied to macroscopic systems. It 
is then reasonable to ask what are the properties of quan- 
tum states of macroscopic systems. In order to fit experi- 
ence, such a quantum state should yield extremely small 
standard deviation for both the position and the mo- 
mentum of the center of mass. On the other hand, due 
to the large number of particles contained in macroscopic 
bodies, they possess an extremely high energy level den- 
sity, or an extremely small average level spacing. Corre- 
spondingly, a macroscopic system M cannot be isolated 
from its environment hence, according to orthodox 
quantum mechanics, it cannot be characterized by a wave 
function, only by a suitable reduced density matrix p, re- 
sulting from the wave function ^ of the composite system 
M + E by taking the partial trace of >< W| over the 
Hilbert space of the environment. Q] 

Nevertheless, it is still possible to assign a wave func- 
tion to the system M. To show this, let us consider a sim- 
plified description of a spin measurement on a spin-half 
particle. Suppose that the S z spin component is mea- 
sured, and initially the particle is in the state | f>, which 
is the eigenstate of the S z operator corresponding to the 
eigenvalue +§• Then the measurement process may be 
symbolically written as | j> |mo >^ | |> >, where 
| mo > stands for the initial state of the measuring de- 
vice, and |m-f > stands for the state of the measuring 
device when it shows the result +6. The horizontal ar- 
row denotes the (unitary) time evolution, satisfying the 



1 The environment itself is a physical system, consisting of 
air molecules, radiation field, other macroscopic systems etc. 



Schrodinger equation. 

Similarly, if the initial state of the particle is | J.>, 
corresponding to the eigenvalue — |, the measurement 
process can be written as | |> |mo > — > | J.> |wj >. 

If the initial state of the particle is the superposition 
"I T> +P\ i> (with \a\ 2 + \(3\ 2 = 1), then, owing to 
the linearity of the Schrodinger equation we may write 
(a| T> +P\ j>) |m >^ a\ f> |m t > +0\ |> \m l >. 

Looking at the final state of the system, one can see 
that the measuring device has no own wave function, 
instead, it can be characterized by the density matrix 
p = |m T > H 2 < m T | + \m i > |/3| 2 < mj. 

Certainly, in such a situation the result of the mea- 
surement is either +§, corresponding to the state |mj >, 
or — §, corresponding to the state |mj, >. As the states 
|m| > and \m± > are orthogonal, they are just the eigen- 
states of p. Extrapolating the conclusion of this simple 
example, it is reasonable to assume that a macroscopic 
system is described by one of the eigenstates of its re- 
duced density matrix p, calculated from the wave func- 
tion of the Universe, and the probability that this par- 
ticular eigenstate occurs is given by the corresponding 
eigenvalue. This has already been proposed by several 
authors H, ||, ffl|. A consistent theory involving this 
idea can be found in Ref. ||]. It is of worth mentioning 
that this latter theory is an explicit counterexample to 
Bell's theorem [|. 

The next question is how the reduced density matrix 
of a macroscopic system looks like if the interaction with 
the environment is taken into account. This is the central 
problem of the phenomenon 'environment induced deco- 
herence' Q, fflb, and as such, has been much studied in 
model systems. A typical model JlO[ contains a number of 
harmonic oscillators playing the role of a heat bath (this 
is the environment) and a macroscopic object coupled to 
it. As a result of the interaction, the reduced density 
matrix of the macroscopic object becomes nearly diag- 
onal in both coordinate and momentum representation, 
so that the width across the diagonal is microscopically 
small, while the distribution along the diagonal (this is 
the probability distribution of the coordinate in coordi- 
nate representation and the probability distribution of 
the momentum in momentum representation) extends to 
a macroscopic regime. Note that (up to the present au- 
thor's best knowledge) none of the decoherence models 
studied so far possesses the property that the eigenstates 
of the reduced density matrix are generically localized. 
The reason may well be that these models are oversim- 
plified. 
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Considering now the mathematical consequences of the 
above idea, one finds that the quantum state of a macro- 
scopic system is a member of an orthonormed set of 
states^. The relevant eigenstates of the reduced den- 
sity matrix (i.e., those which occur with non-negligible 
probability) must correspond to a macroscopic region of 
the classical phase space, therefore, their number is given 
by this phase space volume divided by Planck's constant 
h (in one dimension, in / dimension one must use h' 
instead). All the relevant eigenstates must possess the 
classical property that the standard deviation for both 
the coordinates and the momenta (Sx and Sp, respec- 
tively) is microscopically small. The aim of the present 
paper is to answer the question whether it is possible 
mathematically at all. 

It is natural to ask whether there exists such a com- 
plete, orthonormed basis that the standard deviations Sx, 
Sp for all the basis functions lie below a common bound. 
The answer is known in the one-dimensional case when 
the basis functions are generated from the same function 
by translations along the x and p directions . It turns 
out that completeness, orthonormality and boundedness 
of Sx and Sp cannot be simultaneously fulfilled. 

We made a numerical study to decide whether it is also 
true if the basis functions are not generated by trans- 
lations. The method was the following. Initially one 
dimensional oscillator eigenstates \n > below a given 
energy Eq = (N — l/2)hw were considered. Certainly, 
in this case Sx n cx Sp n cx i/n. Then new basis func- 
tions \n > were constructed by a unitary transformation 

I" >= E^o 1 U nj\j >■ 

The transformation matrix was chosen so as to mini- 
mize the sum of Sx 2 + Sp 2 for all the N basis functions. 
(The units were chosen so that h = 1, u> = 1, m = 1.) As 



Sx 2 —< h\x 2 \n > — < h\x\h > 2 , and Yln=0 < 



n\x \n > 



E^o 1 E-To 1 ELV EGA* < ^ 2 l fc > 



= EjLo 1 < j\% 2 \j > (similarly for p), it was enough 
to maximize the quantity S = E^Lo* ^ H|x|n > 2 
+ < h\p\n > 2 . 

The transformation matrix was built up as a product of 
unitary matrices, each having only one nontrivial 2 by 2 
block. The corresponding indices were randomly chosen. 
For the other indices the matrix was just the unit ma- 
trix. The most general 2 by 2 unitary transformation is 
the transformation matrix of the basic spinors which may 
be parametrized by three angles. For each matrix, these 
angles were chosen at random, and it was subsequently 
considered whether the result of this transformation gave 
rise to a larger sum S or not. If yes, then the transforma- 
tion was kept, if not, it was discarded. The method was 



therefore a simple Monte-Carlo simulation. The proce- 
dure was continued until the quantity S saturated. The 
resulting average deviations (i.e., the sum of 5x 2 + dp 2 
for all the N states divided by N) as a function of var- 
ious N values is plotted in Fig.l. The error of the data 
is less than 10 -2 . The measured values are surprisingly 
well fitted by the logarithmic function 1.0 + 0.6 In (AT). 
This result supports the generalization of the statement 
of Ref. [n], i.e., it shows that no complete, orthonormed 
basis with bounded standard deviations Sx, 5p exists. 




FIG. 1. Diamonds: the average of Sx 2 + Sp 2 over all 
the N states versus N. Dashed line: the logarithmic fit 
1.0 + 0.61og(iV). 

Returning to the original question, we show that the 
logarithmic dependence on N is consistent with the re- 
quirement that the relevant eigenstates (i.e., those with 
significant probability to occur) of the reduced density 
matrix must have small 5x, Sp values^. Indeed, the rele- 
vant eigenstates of p are associated with a large but finite 
portion of the classical phase space. We choose the num- 
ber N to be the ratio of this phase space volume and 
Planck's constant The widths Sx, Sp of the relevant 
eigenstates are proportional to ^/0.5 + 0.3 In (N). This 
quantity remains moderate even for an astronomically 
large N value. As an example, consider the motion of the 
Earth. The linear size in coordinate space is 3 x 10 11 m, 
and that in momentum space is 3.6 x 10 29 kg— . Therefore, 



2 The latter property is due to the hermiticity of the reduced 
density matrix p. 



The irrelevant eigenstates may be identified with the origi- 
nal oscillator eigenstates |n > with n > N. 

4 This is valid in one dimension. In dimension d we divide 
by h d , and the number N is chosen to be the d-th root of the 
ratio, as the basis functions will be constructed as a d-fold 
direct product of the one-dimensional basis. For simplicity 
we shall consider only the one-dimensional case. 
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N = 1.6 x 10 74 and ^0.5 + 0.3 In (Fj = 7.2. Thus the 
widths cix, 5p are only one order of magnitude larger than 
those of the minimal wave packet |0 >. As we see, the log- 
arithmic dependence makes the above assumption about 
the wave functions of macroscopic systems physically ac- 
ceptable. Nevertheless, the problem is left open whether 
these states actually occur in Nature and if so, what the 
origin of their 'localization' property is. Certainly, it 
should be due to the interaction between the macroscopic 
system and its environment (as this produces the mixed 
state described by p), but that feature of the dynamics 
which is responsible for classical behavior is unknown. It 
should be emphasized that the 'nearly diagonal' shape of 
p in coordinate and in momentum representation is insuf- 



'f _ (x+x')' 



with 



ficient. E.g., if p(x,x') oc exp f — 

a i << (72, the density matrix is nearly diagonal in both 
coordinate and momentum representation, but its eigen- 
states are oscillator eigenstates with 8x n oc 5p n oc i/n. 

Here we consider only some consequences if the above 
'logarithmic localization' of the macroscopic quantum 
states actually occurs. 

Let us consider first a harmonic oscillator of unit mass 
and frequency (certainly, these units may be as large as 
1 kg and 1 s _1 ), and calculate the mean deviation of 
the energy in the above localized states. Certainly, the 
average of this quantity depends on the number of the 
states (N), as for larger N the states are broader. The 
dependence of the average square deviation of the energy 
on the number of states is plotted in Fig. 2. 
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FIG. 2. Diamonds: the average of 8E 2 (the squared devi- 
ation of the energy) over all the N states versus N. Dashed 
line: the fit 0.44 iV 125 

Somewhat surprisingly, this dependence is not loga- 
rithmic, but power-like, with the exponent ~ 1.25. This 
means that 



The mean energy is proportional to N (as N is cho- 
sen according to the relevant classical phase space vol- 
ume), therefore, the relative uncertainty of the energy, 



i.e., 5E/E goes to zero as N' 



-0.39 



-0.39 



Certainly, 



in case of a macroscopic oscillator this is unobservable, 
although SE itself can be quite large at the microscopic 
scale fruj. 

In order to understand the behavior (|TJ) , let us consider 
the asymptotic behavior of the above localized states in 
coordinate representation. Beyond \fN they fall off as 
exp (— x 2 /2), due to their construction. For smaller dis- 
tances (down to some xq) they decay like a power. The 
exponent is v w —1.5. This readily implies that 

/oc fvN 
dx \tp(x) | 2 x 2 w const, x / dx x" 1 
- OO J Xq 

= ai + bi In (TV) 

where ai and bi are constants. The behavior in momen- 
tum space is the same. As for the square deviation of 
the energy, we have to consider the expectation value of 
x 4 + p 4 + x 2 p 2 + p 2 i 2 . We may write e.g. 



<x 4 >-- 



dx\ip(x)\ 2 x 4 



const, x / dx x 

J XO 

= a 2 + b 2 A 



One can derive the same linear estimate for < p 4 > as 
well. This implies that 5 2 E grows at least linearly with 
N, which is consistent with the numerically obtained 
N 1 - 25 . 

In order to get further consequences, let us consider a 
macroscopic system being in thermal equilibrium with a 
heat bath at temperature T. What is the density matrix 
of that system at a given instant of time? One would 
guess that it is 



1 



■ exp 



(2) 



SE oc JV' 



0.63 



(1) 



where (3 = an< ^ ^ stands for the canonical par- 

tition function. This expression, however, substitutes 
a long time average, therefore, one cannot expect that 
it is strictly valid at a given instant of time. Accord- 
ing to the previous considerations, it is more realis- 
tic to assume that (||) gives only the long time aver- 
age of the density matrix, which is actually given by 
p = Yin \ L P n > Pn < <Pn\> where \<p n > is narrow in both 
coordinate and in momentum space if p n is not negligibly 
small. As energy eigenstates are generically not local- 
ized, the \ip n >-s are not eigenstates of the Hamiltonian 
H, therefore, p will have nondiagonal matrix elements in 
energy representation. Explicitly, these matrix elements 

are p jk =< Cj|p|Cfc >= E„-P« < Sjl^n >< <Pn\Zk >, 
where |£fc >-s stand for the energy eigenstates. The 
states \ip n > are also localized in energy representation 
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(i.e., for a given n < £,j\tp n > takes on significant values 
only for a narrow range of the energy eigenvalues Ej), 
thus pj.k is a narrow band matrix. 

The existence of nondiagonal matrix elements in en- 
ergy representation can in principle lead to experimen- 
tally observable effects. Suppose one isolates the macro- 
scopic system in question so that the interaction energy 
AEi with the environment becomes much less than the 
effective width SE of the density matrix in energy repre- 
sentation. Under such circumstances consider the linear 
response of the system. For a time t which is much less 
than h/AEj but much larger than h/SE the deviation in 
the expectation value of an operator A is 



5 < A > (t) 



~ / dt'Tr 
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where Hi is the perturbation, and the operators are 
Heisenberg operators calculated with the unperturbed 
Hamiltonian Hq. Choosing A — Hq, the response can 
be cast to the form 

S < H > (t) = ]T ( Pj , k (t) - Pj,k(0)) < &l#i(0)& > , 

where pj,k(t) = exp (-j-(E k - Ej)t) Pj,k(0). If the den- 
sity matrix were strictly diagonal in energy representa- 
tion, the response 5 < Ho > (t) would vanish. Actually, 
as noted above, /Qj,fc(0) is a narrow band matrix, hence 
the Fourier transform of the response vanishes at high 
frequences but becomes significant at low frequences. 
Therefore, the localization property of the eigenstates 
shows up in a suitable linear response which increases 
at low frequencies. It is a challenging question whether 
this behavior may be related to the celebrated 1// noise 
Jlj| . It is also of worth noting that the deviation of the 
density matrix from the canonical expression <Q) implies 
the breakdown of the fluctuation-dissipation theorem. 
In conclusion, we have seen that 

i, in one dimension there is no localized basis in the strict 
sense 

ii, the numerically found 'logarithmically localized' basis 
demonstrates that the idea that macroscopic systems are 
described by the eigenstates of their reduced density ma- 
trix is consistent with experience 

iii, this idea implies that the density matrix of a macro- 
scopic system being in thermal equilibrium is not diag- 
onal in energy representation, that may be observed in 
the spectrum of a suitable linear response function as an 
enhancement at low frequences. 
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